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Abstract 

This is the documentation of the tomographic X-ray data of a wal¬ 
nut made available at http://www.fips.fi/dataset.php, The data 
can be freely used for scientific purposes with appropriate references 
to the data and to this document in arXiv, The data set consists of 
(1) the X-ray sinogram of a single 2D slice of the walnut with three 
different resolutions and (2) the corresponding measurement matri¬ 
ces modeling the linear operation of the X-ray transform. Each of 
these sinograms was obtained from a measured 120-projection fan- 
beam sinogram by down-sampling and taking logarithms. The orig¬ 
inal (measured) sinogram is also provided in its original form and 
resolution. In addition, a larger set of 1200 projections of the same 
walnut was measured and a high-resolution filtered back-projection 
reconstruction was computed from this data; both the sinogram and 
the EBP reconstruction are included in the data set, the latter serving 
as a ground truth reconstruction. 


1 Introduction 

The main idea behind the project was to create real CT measure¬ 
ment data for testing sparse-data tomography algorithms. A walnut 
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shares some of the challenging features with the targets in typical CT 
applications. It contains a dense, layered shell enclosing various struc¬ 
tures with different shapes and contrasts. Like a skull, the shell has 
reflection-symmetry, while the edible part of a walnut contains enough 
non-convexity for testing purposes. 

The CT data in this data set has been used for testing a wavelet- 
based sparsity-promoting reconstruction method and the so-called “S- 
curve method” for choosing the regularization parameter, see [2] or 
|3]. Moreover, the same data was used for testing the total variation 
reconstruction algorithms introduced and studied in [T]. 

2 Contents of the data set 

The data set contains the following MATLABj^data files: 

• Data82.niat, 

• Datal64.mat, 

• Data328.mat, 

• FullSizeSinograms.mat and 

• GroundTruthReconstruction.mat. 

The first three of these files contain CT sinograms and the corre¬ 
sponding measurement matrices with three different resolutions; the 
data in files Data82.mat, Datal64.mat and Data328.mat lead to re¬ 
constructions with resolutions 82 x 82, 164 x 164 and 328 x 328, re¬ 
spectively. The data hie named FullSizeSinograms .mat includes 
the original (measured) sinograms of 120 and 1200 projections, and 
GroundTruthReconstruction.mat contains a high-resolution FBP re¬ 
construction computed from the 1200-projection sinogram. Detailed 
contents of each data hie below. 

Data82.mat contains the following variables: 

1. Sparse matrix A of size 9840 x 6724; measurement matrix. 

2. Matrix m of size 82 x 120; sinogram (120 projections). 

Datal64.mat contains the following variables: 

1. Sparse matrix A of size 19680 x 26896; measurement matrix. 

2. Matrix m of size 164 x 120; sinogram (120 projections). 

Data328.mat contains the following variables: 

^MATLAB is a registered trademark of The Math Works, Inc. 
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1. Sparse matrix A of size 39360 x 107584; measurement matrix. 

2. Matrix m of size 328 x 120; sinogram (120 projections). 

FullSizeSinograms .mat contains the following variables: 

1. Matrix sinograml20 of size 2296 x 120; original (measured) sino¬ 
gram of 120 projections. 

2. Matrix sinograml200 of size 2296 x 1200; original (measured) 
sinogram of 1200 projections. 

GroundTruthReconstruction.mat contains the following variables: 

1. Matrix FBP1200 of size 2296 x 2296; a high-resolution filtered 
back-projection reconstruction computed from the larger sino¬ 
gram of 1200 projections of the walnut (“ground truth”). 

Details on the X-ray measurements are described in Section below. 
The model for the CT problem is 

A* X = m(:), (1) 

where m (:) denotes the standard vector form of matrix m in MATLAB 
and X is the reconstruction in vector form. In other words, the recon¬ 
struction task is to hnd a vector x that (approximately) satisfies 0 
and possibly also meets some additional regularization requirements. 
A demonstration of the use of the data is presented in Section]^ below. 

3 X-ray measurements 

The data in the sinograms are X-ray tomographic (CT) data of a 2D 
cross-section of a walnut measured with custom-built fiCT device Nan¬ 
otom supplied by Phoenix|Xray Systems -|- Services GmbH (Wunstorf, 
Germany). The X-ray detector used was a CMOS flat panel detector 
with 2304x2304 pixels of 50 /rm size (Hamamatsu Photonics, Japan). 
The measurement system is illustrated in Figure and the measure¬ 
ment geometry is shown in Figure 

A set of 120 cone-beam projections with resolution 2304 x 2296 
and angular step of three (3) degrees was measured. Each projection 
image was composed of an average of six 750 ms exposures. The x-ray 
tube acceleration voltage was 80 kV and tube current 200 fiA. See 
Figure for two examples of the resulting projection images. From 
the 2D projection images the middle rows corresponding to the central 
horizontal cross-section of the walnut were taken to form a fan-beam 
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Figure 1: Left: Experimental setup used for collecting tomographic X-ray data of 
a walnut. The detector plane is on the left and the X-ray source on the right in the 
picture. The walnut is attached to a computer-controlled rotator platform. Right: 
Two examples of the resulting 2D projection images. The fan-beam data in the 
sinograms consist of the (down-sampled) central rows of the 2D projection images 
(indicated by arrows in the picture.) 


sinogram of resolution 2296 x 120 (variable sinograinl20 in file Full- 
SizeSinograms.mat). This sinogram was further down-sampled by 
binning and taken logarithms of to obtain the sinograms m in files 
Data82.mat, Datal64.mat and Data328.mat. 

The organization of the pixels in the sinograms and reconstruc¬ 
tions is illustrated in Figure The pixel sizes of the reconstructions 
are 0.513 mm, 0.257 mm and 0.128 mm in the data in Data82.mat, 
Datal64.mat and Data328.mat, respectively. 

In addition, a larger set of 1200 projections of the same walnut 
using the same imaging setup and measurement geometry but with a 
finer angular step of 0.3 degrees was measured (variable sinograml200 
in file FullSizeSinograms .mat). The high-resolution ground truth re¬ 
construction (variable FBP1200 in file GroundTruthReconstruction.mat) 
was computed from this data using filtered back-projection algorithm, 
see Figure]^ 
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Figure 2: Geometry of the measurement setup. Here FOD and FDD denote 
the focus-to-object distance and the focus-to-detector distance, respectively. The 
width of the detector is denoted by W. 
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Figure 3; The organization of the pixels in the sinograms m = [mi, m 2 ,..., mi 207 v]^ 
and reconstructions x = [xi, X 2 , • • •, X];^ 2 ^ in the data in Data82 .mat, Datal64.mat 
and Data328.mat {N = 82,164 or 328). The picture shows the organization for 
the first projection; after that the target takes 3 degree steps counter-clockwise 
(or equivalently the source and detector take 3 degree steps clockwise) and the 
following columns of m are determined in an analogous manner. 
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Figure 4: The high-resolution filtered back-projection reconstruction FBP1200 of 
the walnut computed from 1200 projections. 
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4 Example of using the data 

The following MATLAB code demonstrates how to use the data. 
The code is also provided as the separate MATLAB script file ex¬ 
ample, m and it assumes the data files (or in this case at least the file 
Datal64.mat) are included in the same directory with the script file. 

% Load the measurement matrix and the sinogram from 
7o file Datal64.mat 
load Datal64 A m 

7o Compute a Tikhonov regularized reconstruction using 

7o conjugate gradient algorithm pcg.m 

N = sqrt(size(A,2)); 

alpha = 10; 7o regularization parameter 

fun = @(x) A.’*(A*x)+alpha*x; 

b = A.’*m(:); 

X = pcg(fun,b); 

7o Compute a Tikhonov regularized reconstruction from only 

7o 20 projections 

[mm,nn] = size(m); 

ind = [] ; 

for iii=l:nn/6 

ind = [ind,(1:mm)+(6*iii-6)*mm]; 

end 

m2 = m(:,1:6:end); 

A = A(ind,:); 

alpha = 10; 7o regularization parameter 

fun = @(x) A.’*(A*x)+alpha*x; 

b = A.’*m2(:); 

x2 = pcg(fun,b); 

7o Take a look at the sinograms and the reconstructions 
figure 

subplot(2,2,1) 
imagesc(m) 
colormap gray 
axis square 
axis off 

title(’Sinogram, 120 projections’) 
subplot(2,2,3) 
imagesc(m2) 
colormap gray 



axis square 
axis off 

title(’Sinogram, 20 projections’) 

subplot(2,2,2) 

imagesc(reshape(x,N,N)) 

colormap gray 

axis square 

axis off 

title({’Tikhonov reconstruction,’; ’120 projections’}) 

subplot(2,2,4) 

imagesc(reshape(x2,N,N)) 

colormap gray 

axis square 

axis off 

title({’Tikhonov reconstruction,’; ’20 projections’}) 
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